Quantum noise and mode nonorthogonality 
in nonhermitian PT-symmetric optical resonators 
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"PT-symmetric optical resonators combine absorbing regions with active, amplifying regions. The 
latter are the source of radiation generated via spontaneous and stimulated emission, which embod- 
ies quantum noise and can result in lasing. We calculate the frequency-resolved output radiation 
intensity of such systems and relate it to a suitable measure of excess noise and mode nonorthog- 
onality. The lineshape differs depending on whether the emission lines are isolated (as for weakly 
amplifying, almost hermitian systems) or overlapping (as for the almost degenerate resonances in 
the vicinity of exceptional points associated to spontaneous PT-symmetry breaking). The calcu- 
lations are carried out in the scattering input-output formalism, and are illustrated for a quasi 
one-dimensional resonator set-up. In our derivations we also allow for the more general case of a 
resonator in which the amplifying and absorbing regions are not related by symmetry. 

PACS numbers: 42.50.Nn,03.65.Nk,42.25.Bs,42.55.Ah 
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I. INTRODUCTION 



Recent theoretical and experimental advances in optics 
[l|-[l3j have raised the prospect to realize nonhermitian 
systems with a real spectrum but nonorthogonal wave 
functions. These efforts are based on the concept of VT 
symmetry, originally formulated as a variant of quantum 
mechanics in which potentials can be complex [l3 - [l9j . 
In the optical context such potentials can be realized via 
absorbing and amplifying regions. In a T^T-symmetric 
situation, a discrete unitary operation (such as a reflec- 
tion or inversion) maps the absorbing onto the ampli- 
fying parts, with matching absorption and amplification 
rates. On the level of classical optics, absorption and 
amplification are related by an antiunitary time-reversal 
operation. If the absorption and amplification rates are 
low enough, the spectrum of the composed system is real, 
but beyond a threshold pairs of complex conjugate eigen- 
values appear [l^-^] . This transition (known as sponta- 
neous T'T-symmetry breaking) can be exploited to realize 
a number of exoti c op tical effects, such as unidirectional 
transmission [3!,[8!,'l0], absorption-enhanced transmission 
@ , power oscillations [3, ^] , nonlinear switching H, and 
coexistence of lasing and perfect absorption [l^, ■ At 
the transition point, eigenvalues coalesce, resulting in an 
exceptional point where the two eigenmodes become de- 
generate not only in frequency, but also share the same 
wave function [22h25I |. This singular scenario receives 
considerable attention also for optical systems without 
VT symmetry 

In this paper, we investigate how the unavoidable con- 
sequences of leakage, instability, and quantum noise affect 
the characteristics of realistic T^T-symmetric resonators. 
In combination, we find that these effects offer a window 
to directly access the signatures of nonhermiticity. In 
particular, they allow to detect mode nonorthogonality, 
which discriminates these systems from ordinary hermi- 



tian systems that possess a real spectrum but feature 
mutually orthogonal eigenmodes. 

In realistic devices, additional losses arise due to leak- 
age, as radiation needs to be coupled out of the system. 
Even though these losses break exact T'T-symmetry, sig- 
natures of the associated peculiar spectral characteris- 
tics are still present in the complex resonance frequencies 
of the open system. The coalescence of eigenvalues at 
the spontaneous T^T-symmetry breaking transition then 
translates to situations where two resonance frequencies 
approach each other very closely in the complex plane. 

In actively amplifying optical systems, the appearance 
of real eigenfrequencies indicates an instability, the on- 
set of lasing. The consequences for T'T-symmetric sys- 
tems have been explored only since very recently. In 
these systems, the lasing threshold is either reached in 
the limit of the closed system (if the spectrum in this 
limit is real) or at finite leakiness (if VT symmetry 
in the closed system is spontaneously broken, i.e., beyond 
an exceptional point) [H, In both cases, the system 
is in practice stabilized by saturation in the amplifying 
parts, thereby assuring that the output intensity remains 
finite. This saturation provides a physical mechanism 
that breaks the balance of amplification and absorption 
required for VT symmetry. 

An ordinary laser emits coherent radiation with a nar- 
row emission line that can be well approximated by a 
Lorentzian. According to general laser theory [s^, [33| . 
the width Aw of the Lorentzian arises due to noise, of 
which a certain amount, quantum noise, is an unavoid- 
able consequence of the quantum nature of microscopic 
emission events. Investigations on purely amplifying (not 
T^T-symmetric) systems established a direct link of non- 
hermiticity and an enhanced line broadening (known as 
excess noise), which are both captured by a measure of 
mode nonorthogonality, the Petermann factor K (34l439j . 
At an exceptional point, K diverges because of the coales- 
cence of resonance wavefunctions [23 - l25j . PT-symmetric 
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systems offer an ideal venue where the consequences for 
the radiated intensity in this singular case can be ex- 
plored. More generally, one should expect for such sys- 
tems that the excess noise provides a probe of the level 
of nonhermiticity also away from an exceptional point. 

The preceding observations capture our principal mo- 
tivation for this work. It is the purpose of this paper to 
formulate a theory of the quantum noise and radiation 
of leaky T^T-symmetric optical systems in the full range 
of situations far below, near, and beyond the reconfig- 
uration of the spectrum at an exceptional point, up to 
the point where the lasing threshold is reached. This re- 
quires a quantum optical treatment, which we base on 
the scattering input-output formalism |40l442i . as previ- 
ously applied to purely amplifying systems |38l , |39| . Tak- 
ing the absorbing parts of the system into account, we 
establish general relations for the output intensity as for 
the previously studied case of a homogeneously ampli- 
fying resonator, but find that this involves a nontrivial 
combination of aspects from mode nonorthogonality and 
excess noise. As one approaches an exceptional point, the 
partial intensities of the two near-degenerate resonances 
still diverge, but the combined amplitude remains finite, 
which signals a change in the lineshape from a Lorentzian 
to a squared Lorentzian (as also observed at exceptional 
points in passive scattering theory (43l - |45| ). 

We illustrate these general results on quantum noise 
for a specific quasi one-dimensional P7~-symmetric res- 
onator, which displays the generic spectral properties 
of previously studied resonators [H, and offers ad- 
ditional control via a variable leakage to the exterior. 
This application extends the investigation in Ref. pT| . 
which used the input-output formalism to study a VT- 
symmetric resonator with well isolated resonances and 
did not address the relation to mode nonorthogonal- 
ity. Throughout our derivations, we also present gen- 
eral expressions that apply to systems with amplify- 
ing and absorbing regions, even when these are not re- 
lated hy VT symmetry (as recently investigated, e.g., in 
Refs. B Hi], El 113) ■ 

This work is organized as follows. Section |lT] reviews 
the spectral features of closed and open T^T-symmetric 
systems, as well as the signatures of excess quantum noise 
for conventional, amplifying resonators, with the discus- 
sion based in both cases on the common framework of 
(first- and second-quantized) scattering theory. In Sec- 
tion |lll] we adapt the scattering input-output formalism 
to resonators with absorbing parts and derive general ex- 
pressions for the output radiation intensity. In the cen- 
tral Section ITVl we analyze this radiation near resonance 
and establish the link to mode nonorthogonality. Sec- 
tion |V] sees our general results applied to a specific VT- 
symmetric resonator set-up. We first study the classical 
wave problem and determine the resonance frequencies 
and exceptional points. We then analyze the output ra- 
diation in the vicinity of these frequencies and verify the 
sensitivity to the nonorthogonality of modes, as well as 
the emergence of a squared Lorentzian at the exceptional 



points. Section IVll contains our conclusions. 

II. BACKGROUND AND OPEN ISSUES 

In this section we briefly review the ('first-quantized', 
classical- wave) scattering approach to the determination 
of eigenfrequencies and resonances in closed and open 
nonhermitian 'PT-symmetric systems, as well as the ap- 
plication of the ('second-quantized', quantum-optical) 
scattering input-output formalism to the problem of ex- 
cess noise in purely amplifying systems. Sections IIIII and 
llVl then extend the latter to partially absorbing systems, 
including those with VT symmetry. Throughout, we 
identify uj with the energy of photons (effectively setting 
h= 1). 

A. Scattering quantization and spectral properties 
of "PT-symmetric systems 

While scattering appears most naturally in the context 
of transport through open systems, the scattering frame- 
work can also serve as a convenient technique to deter- 
mine spectral properties, such as the resonance eigenfre- 
quencies of these systeins, or even the bound-state spec- 
trum of closed systems [48|, fio^ . In scattering theory the 
eigenfrequencies are determined efficiently using a small 
number of basis functions, which are the on-shell (fixed 
energy or frequency) solutions in various, artificially sep- 
arated (and therefore open) parts of the system, which 
then are coupled together using their scattering matri- 
ces. The quantization condition follows from the con- 
sistency requirement of the matching conditions. In the 
present work, the use of scattering theory is further mo- 
tivated by the fact that it also can serve as a framework 
for introducing quantum noise. Indeed, we will find that 
the steps in the derivation of the quantization condition 
(given here) and in the implementation of quantum noise 
in composed systems (given in Sec. IIIip closely resemble 
each other. In this first part of background material, we 
therefore review the general ideas of spectral analysis in 
scattering theory, and also describe the consequences for 
systems with VT symmetry [ll|, [l^ (for other aspects of 
PT-symmetric scattering theory see Refs. |50l - l55j ). 

1. Scattenng matrix 

We start with the definition of the relevant scatter- 
ing matrix. Consider a classical wave equation, e.g., the 
Helmholtz equation 

A?/'(r) + c-^uj^n'^{r)ij{r) = (1) 

for TM polarized light in a two-dimensional dielectric 
resonator, where scattering, absorption, and amplifica- 
tion enter via the dielectric index n(r) (here c is the 
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speed of light in vacuum). Assuming Ren > 0, we have 
Im n < in amphfying regions, while in absorbing regions 
Imn > 0. The wave equation is then solved for given 
amplitudes a™ of incident propagating modes, which de- 
livers a linear relation 



(2) 



for the amplitudes a°"* of outgoing modes. Here S{uj) is 
the scattering matrix, whose dimensions depend on the 
number of incoming and outcoming propagating modes 
at the given frequency uj. For a system where modes are 
coupled in from a left or right entrance, the scattering 
amplitudes can be collected into reflection blocks r and 
r' , as well as transmission blocks t and t' (where the 
prime discriminates whether the incident radiation comes 
from the left or right, respectively), such that 



S 



t' 



(3) 



In the special case that amplification and absorption 
are absent (so that the refractive index n is real), and 
if the frequency is real as well, the scattering matrix is 
unitary, S\uj)S{lo) = 1. This relation embodies the con- 
servation of particle flux. For a complex refractive index, 
however, this conservation law is in general violated. 



2. Spectral properties of closed systems 

Consider a closed resonator that is split (if only artifi- 
cially) into two parts L and R ('left' and 'right'), which 
are described by scattering matrices identical to the re- 
flection blocks r^(i^) and ^^ (a;), respectively [see the cen- 
tral region in Fig.[TJa), ignoring for the moment any leak- 
age from the system]. At the interface, amplitudes of 
modes travelling 'to the right' (i.e. from L to R) are col- 
lected into a vector a^, while those for modes travelling 
into the opposite direction are collected into a vector ag . 



The matching conditions 



are consistent if 



det [r'^{uj)rR{uj) - 1] = 0. 



(4) 



(5) 



For a hcrmitian system, where the scattering matrices are 
unitary, this condition only admits real frequencies. In 
presence of absorption or amplification, however, where 
the scattering matrices are nonunitary, the eigenfrequen- 
cies generally are complex. In both cases, these solutions 
are identical to the eigenvalues of the (possibly nonhcr- 
mitian) underlying Hamiltonian. 

The interest in VT symmetry for nonhermitian sys- 
tems arises because this provides a mechanism where at 
least part of the spectrum can still be real. In the scat- 
tering formalism, this is embodied by the relation jTll ITsj 



rn{u:)=Vrr'^{u:) = [{r'^{L.*)r 



(6) 
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FIG. 1. (Color online) Illustration of the wave-classical and 
quantum-optical treatment of PT-symmetric systems, which 
are composed of an absorbing and an amplifying region. Panel 
(a) defines amplitudes used in the scattering approach to de- 
termine the resonance frequencies in the wave-classical limit 
Panel (b) defines the field operators for the 



(see Sec. Ill A|) 

quantum-optical input-output formalism of Sec. IIIII In panel 
(c) we sketch the quasi one-dimensional model resonator stud- 
ied in Sec. [V] This system is terminated by semitransparent 
mirrors with transmission probability T. In the interior, two 
regions of equal length L/2 are equipped with a refractive in- 
dex riL (left region) and ur (right region), which is constant 
throughout each respective region. VT symmetry is realized 
for riR = n|,. We write ul ~ no{l + ia), where a controls the 
degree of nonhermiticity. 



between the scattering matrices, where the labels L and 
R now refer to the symmetry-related subsystems. The 
quantization condition can then be written as 



det [K (..)}* -r^(u;*)]=0. 



(7) 



On the real frequency axis, this reduces to the condition 

El 



detlmr^(a;) = 0, 



(8) 



which can be generically fulfilled by varying a single real 
parameter (i.e., w), as it involves the determinant of a 
manifestly real matrix. However, the quantization con- 
dition can also have complex solutions, which then occur 
in complex-conjugated pairs. 

The transition between both situations involves pairs 
of frequencies that coalesce on the real axis, and then 
move into the complex plane, where they remain related 
by complex conjugation |14| - [2l| . Typically, this transi- 
tion is driven by increasing the nonhermiticity (the trend 
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does not need to be strict; sometimes pairs of resonance 
frequencies become real when the nonhermiticity is in- 
creased). However, the transition also depends on the 
coupling of the absorbing and amplifying regions, and in- 
deed can be induced by reducing this coupling |53i,i56|,i57i] . 
(In the trivial limit where the two regions are decoupled, 
they both possess a fully complex spectrum.) 

While VT symmetry leaves wavefunctions of real 
eigenvalues invariant, it interchanges those of complex- 
conjugated pairs, which are therefore not 'PT-symmetric 
when taken each on their own [l3, [13 ■ This entails that 
at the point of degeneracy the two wavefunctions have to 
collapse, thereby resulting in an exceptional point, which 
is the generic degeneracy scenario in nonhermitian sys- 
tems [23 - l25| . For the case of the Helmholtz equation, 
this behavior can be quantified on the basis of the bi- 
orthogonality relation 



n ipiil;2dr = 0, 



(9) 



which holds for any two resonance modes ■01 j "02, even 
if the refractive index is complex. As one approaches 
the exceptional point, '02 — = "0 is shared between 
the two degenerate eigenvalues, and the wave function 
becomes self-orthogonal. 



n ip dr — (exceptional point). 



(10) 



3. Spectral properties of open systems 

For a leaky (geometrically open) system, such as a res- 
onator confined by semitransparent mirrors, the eigenfre- 
quencies are generally shifted into the complex plane, cor- 
responding to resonance frequencies of quasibound states 
with decay rate — 2Ima;. These quasibound states fulfill 
the wave equation with purely outgoing boundary condi- 
tions. For a passive system, their decay rates are positive. 
In the presence of amplification counteracting the leak- 
age, individual resonances can cross from the lower to the 
upper half of the complex plane. This signifies an insta- 
bility, which in optics corresponds to the laser threshold 

The complex resonance spectrum can again be ob- 
tained from scattering theory Q . Including the leakage 
channels to the outside, the scattering matrices of the 
left and right parts assume a natural block structure 



Sl 



tL r' 



Sr- 



tR r' 



(11) 



which now also contains scattering amplitudes related to 
reflection and transmission from and to the exterior re- 
gions [i.e., from the left and right entrances into these 
subsystems; see Fig.[lja)]. The internal matching condi- 
tions still only involve the reflection blocks of the left 
system and tr of the right system, and thus remain of the 



form (|4]). Hence, the quantization condition is still given 
by Eq. ([S])- However, for the open system this condi- 
tion typically admits only complex solutions, even if the 
underlying Hamiltonian is hcrmitian or T^T-symmetric. 

From the perspective of the open system, the signifl- 
cance of Eq. ([S]) becomes clear when considering the com- 
posed scattering matrix S = Sl° Sr of the total system, 
which is formed according to the matrix composition rule 



tr t'j^ 

tR r' 



rL+t'^j^^rRtL 



:tL 



L l — rBX'i^ ^ 

-tR-rA — r'rt' 



(12) 



These expressions contain the resonant denominators 1 — 
r'j^rR, and thus diverge when condition ([5]) is met. 

In the 7-*T-symmetric case, the scattering matrices in 
Eq. pT|) are furthermore related according to [ll| 

Sr{u) = VTSLiio) - a,[52(c^*)]"V, 

1 r^'-i+. 1 1 



1,1 1 



' T t j — 



1 



,(13) 



where in the flrst line is a Pauli matrix, and the second 
line follows from block inversion formulas. The resonance 
condition can then be written in the form 

det(ri^(..)-[r^(c.*)-t4c.*)r^\..*)ti^(c.*)]*) =0. (14) 

The quantization condition in the closed system is 
recovered for — t'j^ — 0. 

When the system is almost closed, the eigenfrequen- 
cies are only slightly shifted downwards into the com- 
plex plane. The transition of the spectrum from real 
to complex is then replaced by a spectral rearrangement 
in which complex resonance frequencies approach each 
other very closely by moving together roughly parallel 
to the real axis, and then depart from each other into 
the direction of the imaginary axis [1, [H, [I^l- When 
two complex frequencies meet in an exact degeneracy, 
one again reaches an exceptional point, where the wave- 
functions of two resonance modes collapse and become 
self-orthogonal, in accordance to Eq. PU)). 



B. Excess noise for purely amplifying resonators 

A practical method to probe the spectrum of a res- 
onator is to fill it with an active, amplifying medium 
and observe the ensuing emitted radiation intensity. Far 
above the lasing threshold, one then observes coherent 
radiation with a narrow Lorentzian emission line 13211 



^(^) — ^totaiT — 



wo)2-fAwV4' 



(15) 



where the line width Ao; is dictated by noise. For the 
reference point of an almost closed purely amplifying res- 
onator (the prototypical good- cavity laser), the quantum 
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limit of the line width is given by the Schawlow-Townes 
relation [SJ 



Auj = rV2/i 



total 1 



(16) 



where /total = / /(w) dw, while F is the cold-cavity de- 
cay rate in the passive resonator [corresponding to a real 
refractive index in Eq. ([Ij]. This noise is due to spon- 
taneous emission events, which are incoherent and result 
in phase fluctuations, while amplitude fluctuations are 
suppressed by the feedback from the medium [s^ [36| . In 
the linear regime just below the lasing threshold, on the 
other hand, where the emitted radiation is incoherent, 
amplitude and phase fluctuations are both present and 
equal each other, so that the Schawlow-Townes relation 
reads 



Au = //total- 



(17) 



Both versions of the Schawlow-Townes relation require 
absence of absorption and assume orthogonal resonator 
modes, so that crosstalk between photons emitted into 
different modes can be ignored. In order to identify cor- 
rections due to the violation of these assumptions, it is 
advantageous to concentrate on the linear regime just 
below threshold, where the emitted intensity can be cal- 
culated conveniently in the scattering input-output for- 
malism [40l - l42j . Results can be translated to the lasing 
regime by the simple rule that there the quantum-limited 
line width is still reduced by a factor of two 36i] . In the 
remainder of the background section we collect results for 
the previously studied case of purely amplifying systems 
[34l439| . focussing only on some key points as the details 
of the derivation are encompassed by our more general 
considerations in Sees. Illll and IIVI 

The formalism starts with the classical wave equation, 
e.g., the Helmholtz equation ([1]). For a passive system 
(with a real refractive index) the scattering matrix S{uj) 
is unitary, SS'^ = 1, but in presence of amplification 
this condition no longer holds. A physical consequence 
is the appearance of quantum noise. When amplitudes 
are promoted to annihilation operators a'"'°"*, Eq. ([2]) 
with nonunitary S is not compatible with the bosonic 
commutation relations for both sets of operators. This 
problem can be fixed by introducing auxiliary bosonic 
operators b, such that 



S'(cj)a'" + Q(cj)b'f. 



(18) 



The commutation relations for a'" and a°"* are now com- 
patible if 55^ — QQ^ = 1. This constraint has the status 
of a fluctuation-dissipation theorem. It admits solutions 
when Q is equipped with at least as many columns as 
rows, but does not fix Q completely. However, subject 
to reasonable physical assumptions, the constraint is suf- 
ficient to calculate the output intensity of the system in 
terms of its classical wave scattering properties. 

For illustration we consider the case of total population 
inversion ((b^b) = 0) and no incoming radiation. The 



frequency-resolved output intensity is then of the form 

M 

= i-(a™'ta°"t) ^ J_trQQt = -Itr (^^t - 1). 

ZTT ZTT ZTT 

(19) 

The final expression no longer contains the matrix Q and 
therefore can be calculated purely based on the classi- 
cal wave problem. In the vicinity of resonant emission 
frequencies, where the scattering matrix diverges [see 
Eqs. (O and ([T^]. the radiated intensity is large. Lin- 
earizing in the deviation from the resonance condition, 
one generally finds a Lorentzian line shape 



K 



27r (a;-a;o)2 + Aw74' 



(20) 



which features an extra factor K so that the Schawlow- 
Townes relation Eq. p7)) is replaced by 



AlO = /fFV/total. 



(21) 



The factor K is known as the Petermann factor 13411 . and 
can often be related to mode nonorthogonality |35l439l |. 
For the Helmholtz equation ([1]), this factor takes the ex- 
plicit form [3^ 



K = 



j \'iP\HwL{n^) dv 



J ■i/;2Im(n2) dr 



(22) 



which holds also in presence of inhomogeneity in the re- 
fractive index and the gain. For a homogeneously am- 
plifying resonator, this reduces to Siegman's original ex- 
pression [111 



K = 



|/V2dr|2 



(23) 



The connection of the excess noise to mode nonorthog- 
onality becomes apparent when one considers the combi- 
nation (see also Sec. ITV)) [39i] 



J |V'|^Im(n2) dr 



J ifj'^n'^ dr 



(24) 



which appears in the numerator of Eq. ([201) . The self- 
orthogonality condition (jlOp implies that this combina- 
tion diverges at an exceptional point. In general, the inte- 
gral in the denominator of Eq. (p4)) can be interpreted as 
a suitably weighted overlap integral of the left eigenfunc- 
tion "0 and the right eigenfunction -0*, which represents 
the scalar product between these two types of eigenfunc- 
tion (the stated simple relation of left and right eigen- 
functions holds because the potential in the Helmholtz 
equation is scalar, but breaks down in the presence of 
vector potentials ^^). Furthermore, expressions ([^ 
and ((23)) both imply /sT > 1, where K = \ \s only ob- 
tained in the limit of a passive, closed system, in which 
the wave function is real. In more physical terms, on 
the other hand, expression (P^ signifies that the excess 
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noise is generated in the amplifying parts of the system, 
since the passive parts (with real refractive index) do not 
contribute to the integrals. 

Clearly, this duality of interpretations of the Peter- 
mann factor in terms of mode nonorthogonality and ex- 
cess quantum noise breaks down in presence of absorbing 
parts of the system. The mode nonorthogonality is then 
also induced by the absorbing regions with a positive 
imaginary part of the refractive index, but these on their 
own do not create any radiation. One of the main goals of 
the present work is to identify the role of the Petermann 
factor in the presence of such absorbing elements. We fo- 
cus on PT-symmetric systems, as their peculiar spectral 
properties imply that one can also easily steer close to an 
exceptional point, and examine how the total output in- 
tensity is regularized despite a divergent Petermann fac- 
tor. In the remainder of this paper, we therefore first 
set up a theory of quantum noise in simultaneous pres- 
ence of absorbing and amplifying regions, and then iden- 
tify signatures of nonhermiticity, mode-nonorthogonality, 
and exceptional points in the output radiation intensity, 
followed by an illustration for a specific PT-symmetric 
system. 



III. INPUT-OUTPUT FORMALISM IN THE 
PRESENCE OF ABSORBING REGIONS 

In this section we extend the scattering input-output 
formalism described in Sec. Ill Bl to systems which com- 
bine amplifying and absorbing regions. 

In simultaneous presence of gain by amplification and 
loss by absorption, the input-output relations ([T5)) mod- 
ify into [12 

a°-* = 5(c.)a'" + Q(,) (w)b(,) + Q(,) {u;)hl^ , (25) 

where (/) refers to the absorbing (lossy) regions and 
{g) refers to the amplifying regions (with gain). The 
commutation relations now deliver the constraint SS'' + 
Q{i)QIi) ~ ^{g)^lg) ^ ^' ■^hich no longer uniquely re- 
lates the coupling strengths (encoded in the combinations 
QQ^) to the deviation of the scattering matrix from uni- 
tarity. 

In order to circumvent this problem, we assume that 
the absorbing and amplifying regions are spatially sep- 
arated. For T'T-symmetric systems, this assumption is 
rather natural, as the V symmetry is usually of a geo- 
metric nature, such as a refiection a; — > — x about a plane 
perpendicular to the x-axis. To be specific, we assign 
absorption to the left part of the system and amplifica- 
tion to the right part of the system [see Fig. [Ijb)]. This 
strategy, briefiy sketched and applied to a special case in 
Ref. works whenever a similar separation into dif- 
ferent regions is present due to the physical composition 
of the system, even in absence of symmetries. We there- 
fore consider this slightly more general case, i.e., we first 
only assume that the absorbing and amplifying parts can 



be separated, and only then impose symmetry relations 
between them. 

The classical wave scattering from the right and left 
parts is described by the scattering matrices Sl and Sr, 
whose natural block structure given in Eq. PT|) . The 
corresponding input-output relations read 

where the various vectors of operators are defined in Fig. 
[IJb). Based on their commutation relations, the coupling 
strengths to the medium, 

(gf )(QLq2) = 1-^l4, (27a) 

(§f )(QL0fl)=^«4 -l' (27b) 

are now again related to the scattering properties of the 
purely absorbing or amplifying parts of the system. As in 
the derivation of Eq. (IT9|) , the relations ^7} are sufficient 
to calculate the output radiation of the system in terms 
of the classical scattering properties encoded in Sl and 
Sr. 

Starting from the relations ([25]) . we algebraically elim- 
inate the auxiliary operators and at the interface 
between the two regions, 

= ^-^[tR^R + TRtL^'l + Qflb|j + r^g'^^bi], 

(28a) 

a? = , _ \ [tL^'l + ^L^a'i? + Q'l^l + r^g^bt^]. 

(28b) 

The annihilation operators for outgoing radiation then 
follow by substituting these expressions into 

ar = TLS^l -f g^bi -I- i'^a^, (29a) 
a??"* = rRk'R + Q'r^r + tRS^. (29b) 

We are interested in the radiation intensity which orig- 
inates from the system, i.e., in absence of external incom- 
ing radiation, and therefore demand 

(C^atTL) = &^:r) = (30) 

for all incoming modes n (expectation values of all 
crossterms between different incoming modes also van- 
ish). Within the medium, we assume mode- independent 
expectation values 

(bl^b„L) = /l, (b„^bl^) = fR, (31) 

which can be associated to the excited-state occupa- 
tions gL = /l/(1 + 2/i) in the absorbing and gR = 
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/i?,/ — 1) in the population-inverted amplifying parts 
of the system, respectively. 

Under these conditions, the frequency-dependent out- 
put intensity of the composed system, resolved depending 
on whether it eventually emerges from the absorbing or 
the amplifying part of the system, takes the form 
1 



-outt 

Zn 



= ^tr 



2tt 

h 
27r' 



1 - rn/r 



tQr 



tr 



Ql + t'L 



1 



1 - TBTj^ 
1 



t^rQ'l 



1 



trTl 



(32) 



Ir{^) 



1 

Jl, 
— tr 
27r 



-outt 



VR 



-Q'l 



it 



tR- 



t'lTr 



Q'l 



+ 



Ir 
27r" 



tr 



Q'r. + ^R 



1 



'lQr 



1 - r'^rn^ 
^ t'lQr 



tR- t'^Qr , (33) 

1 - ^L^R 

where all combinations of Q matrices are completely de- 
termined via the relations p7p . These expressions imply 
the appearance of narrow emission lines for frequencies 
close to resonance, which occurs whenever the denomina- 
tor 1 — r'j^rR vanishes, in accordance to the classical- wave 
quantization condition 

For a quasi one-dimensional resonator, the transmis- 
sion and reflection matrices reduce to complex numbers. 
Equations ([5^ and ([55]) then yield the following more 
compact expressions, 

1 \t'L\H\rR\' + \t'R\'-l) 



II 



2tt 
Ztt 



II 



rRr'L? 



\tRr'L?{\rR\ 



\t'R\' 



1) 



1 



\tR\' 



-2Re- 



^rt'l irRt*j 



'-R' R 



(34) 



(35) 



1 - r^TR 

where we now made explicit use of Eq. I\21\ and assumed 
idealized conditions with fj^ = (ground state popula- 
tion in the absorbing parts) and fR = l (total population 
inversion in the amplifying parts). For 7^7~-symmetric 
resonators the amplitudes are again related by Eq. ()13|) , 
and the resonance condition is given by Eq. (I14p . 



IV. NEAR-RESONANT RADIATION 
INTENSITY AND THE ROLE OF THE 
PETERMANN FACTOR 



Starting from Eqs. and (|35l) for the emitted in- 
tensity from a quasi one-dimensional resonator, we now 



evaluate the total intensity I = II + Ir close to reso- 
nance, cj « wq = Rewm, where the complex resonance 
frequency uim fulfills r'j^{uJm)rR{LUm) = 1 according to 
Eq. ([SJ. Keeping only resonant terms we have 

w 1 \tR\H\t'L? + \tRr'L\')\rR.\' + \tRf^l ,^^, 

n^) - TT- n — 17—12 rm i.*^"-' 



27r 



\l-r'^rR\ 



\tR? 



where we grouped the terms for later convenience. Our 
goal is to relate this expression to properties of the 
resonance wave function which fulfills the one- 

dimensional Helmholtz equation 



(37) 



(with km = ^m/c) subject to purely outgoing boundary 
conditions. The derivation proceeds by a number of tech- 
nical steps which lead to a compact final result, Eq. P5)) 
below. 



A. Near-resonant radiation intensity 

For convenience, we normalize the resonance wave 
function such that in the free space (n = 1) to the right 

— 1/2 

of the resonator -i/),,! (x) — km exp{ikmx). Furthermore, 
we also insert a zero-width layer with n = 1 at the in- 
terface between the right and the left region, where the 
wave function then takes the form 



1pm = tjj^ exp{ikmx)kj^^^ + rRtj^ eyip{-ikmx)kj-^'^ . 

(38) 

Here we expressed the amplitudes by the elements of 
the scattering matrix Sr^ Eq. Furthermore, in 

terms of elements of the scattering matrix Sl, the out- 
going component in the free space to the left of the res- 
onator is given by %lJm{x) = akm^"^ ex-p{—ikmx), where 
a = t'l^/tRr'i^ — TRt'j^/tR (the latter equality follows from 
the resonance quantization condition). 

We now evaluate a number of integrals using Eq. (|37p . 
integration by parts, the condition Rewm ^ |Imw,„|, 
and boundary terms matching the stipulated outgoing 
wave amplitudes. Below, C denotes the left part of the 
resonator, and TZ denotes the right part of the resonator. 
We then have 



lu,{kln')\^m? 



n 



n 



Im(V'mV'm) = Im(V'mV'm)loTC 



\rR? + \tR\^ - 1 

\tR? 



(39) 



Analogously, upon extending the integral over the 
whole resonator 



2 2,.,, ,2_\t'L? + \tRr'L? 



Im/Cj^n \'ip. 



n+c 



\tRr'L? 



(40) 



Finally, we employ a similar integral to investigate the 
behavior of the wave function close to resonance, uj = 
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fee « Wo- Due to the detuning, the wave function then 
also possesses incoming components, but as long as the 
detuning is small the outgoing components remain ap- 
proximately unchanged. Therefore, we can assume that 

— 1/2 

to the left of the resonator, ip{x) ^ km [a e'xp{—ikmx)-\- 
P e.yi\)(ikmx)] with a as given before, while to the right 



-1/2 



exp(ifcma;) + 7 exp(— ifcrnS;)]. 



The incoming components of the wave function can 
now be extracted via the following sequence of steps 
(which starts by linearizing fc^ around resonance): 



2 ; 2 

n+c 

2 7„2 \„2„;,2 



2 kfii (^k k' 

{e - ktJn% 



TZ+C 



TZ+C 



TZ+C 

2i(a/3-f 7)- 



(41) 



Furthermore, in terms of the scattering matrix of the 
whole system we have rj3 + t'^ = a and tfi -\- r'7 = 1, 
which determines the coefficients /3 and 7. We can then 
employ the scattering matrix composition rules (|T2|) to 
express these coefficients in terms of the elements oi Sl 
and Sr. Close to resonance, this reduces to 

(a/3 + 7) ^ -l/r « -(1 - r^rfl)/(r^4). (42) 

Based on expressions (HUl) . and (|1T|) . the near- 

resonant radiation intensity p6p can now be rewritten 
as 



/(w) 



1 /■Rlm^m'^^IV'mP/T^+^Iniw^n^lV;, 



(w- a;o)2 Aa;2/4' 



(43) 



This is the central general result in this work. It ap- 
plies to systems with amplifying and absorbing regions, 
including T^T-symmetric and purely amplifying systems, 
as is discussed in detail in the next subsection. The last 
factor is a Lorentzian of width Aw — |2Imw„i|, cen- 
tered at Wo — ReWm- The remaining combination of 
integrals resembles the expression (j24l) for a purely am- 
plifying system. However, one of the integrals in the 
numerator extends over the whole system (where her- 
miticity is broken), while the other is restricted to the 
(amplifying) right part of the resonator. Therefore, this 
expression presents a mixture of the dual interpretations 
of the conventional Petermann factor as a measure of 
mode nonorthogonality and excess noise. Furthermore, 
as in Eq. (j24l) the denominator involves the appropriate 
overlap integral, which diverges at an exceptional point 
as a consequence of the self-orthogonality relation ([TT 



B. Purely amplifying versus "PT-symmetric systems 

In order to get a hold of the general features encoded in 
the near- resonant intensity (1431) , we first describe how one 
recovers from this expression the Petermann factor (j22p 
for a purely amplifying resonator. In such a system, the 
breaking of hermiticity via the gain (encoded in Im n < 
0) is constrained because this systematically shifts the 
resonances upwards in the complex plane. Therefore, the 
system is driven to the lasing threshold, which is reached 
when a resonance approaches the real axis. For highly 
excited modes (with large wo), this happens very quickly, 
for |Imn| ~ (r/wo)Ren <C Ren, where F is the cold- 

jerturbation theory to 

m 



cavity rate. One then can use 
relate the required gain to F via 



F 

Wo 



/Im (w^n2)V'2 



(purely amplifying res- 
onator at threshold). 

(44) 

Here, the refractive index is no longer arbitrary but has to 
be chosen such that the resonator is at threshold. With 
the help of this relation, Eq. (gS]) turns into Eq. (pOj) . 
recovering the Petermann factor K as defined in ([55]) . 

For systems that are not purely amplifying, on the 
other hand, the breaking of hermiticity does not cause 
a systematic shift of resonances, as the absorbing regions 
counteract the effect of the amplifying regions. In com- 
parison to the Lorentzian (|20l) . one could then define a 



Wn 



generalized Petermann factor as 



K 



Itz Im (w2„n2)|V'„^P Im (w^?^^)|V'. 



Y2- 



(45) 

The same perturbative treatment that leads to Eq. pi)) 
entails more generally that the cold-cavity decay rate F 
is related to the (measurable) line width Aw = — 2Imw 
according to 



F 

Wo 



Aw 

Wo 



/Im(w2^n2)V'2 



(absorbing & amplify- 
ing resonator, not nec- 
essarily at threshold). 

(46) 

[Equation (H^ follows by demanding Aw ^ F close to 
threshold.] 

For the specific case of a PT-symmetric resonator, 
the validity of this perturbative treatment is limited to 
the regime before the exceptional point, as such a point 
induces degeneracy, and the self-orthogonality property 
([TU| implies that the denominator in Eq. (|46p diverges. 
Before one reaches the exceptional point, however, the 
individual wave functions are T^T-symmetric, so that the 
integral in the numerator of Eq. (j46p vanishes. Therefore, 
one can safely approximate 



F = Aw (T'T-symmetric resonators). 



(47) 



which in practice should hold up to very close to the ex- 
ceptional point. Beyond the exceptional point, it should 
initially be reasonable to approximate F by the average 
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width the two involved, overlapping resonances. Typ- 
ically (and as we will confirm below), the regime far 
beyond the exceptional point is physically inaccessible 
as one quickly reaches the laser threshold [Ts"]. Based 
on these observations, Eq. (|43p can be directly applied 
to specific T'T-symmetric systems, which we illustrate 
in the next section for the example of a quasi one- 
dimensional resonator setup. 



V. APPLICATION TO A PT-SYMMETRIC 
RESONATOR 

A. Model system set-up 

We now consider a specific T'T-symmetric resonator 
[depicted in Fig. HJc)], which is made of two regions 
of equal length L/2 (i.e., the total resonator length is 
L). In the left (absorbing) part of the system, the com- 
plex refractive index = + {'"^i > 0)i while 
UR = n*j^ = no — ini in the right (amplifying) part of 
the system. The resonator is terminated by two identi- 
cal semitransparent mirrors of transmission probability 
T. This resonator can be interpreted as an open ver- 
sion of the system studied in Ref . flal . It is also similar 
to the resonator studied in Ref. [llj, but features some 
backscattering at the interface between the regions due 
to the step in Imn. Furthermore, apart from the ad- 
ditional mirrors, it is like the resonator studied in Ref. 
(l3i] . These slight modifications are motivated by specific 
requirements for our investigation. The system must be 
open to study the output radiation, the backscattering 
facilitates the appearance of exceptional points, and the 
mirrors allow to study the limit T — >• of an almost closed 
resonator. Furthermore, the degree of nonhermiticity can 
be controlled by changing nj. 



B. Resonance frequencies and exceptional points 

In order to obtain the resonance frequencies we ap- 
ply the scattering quantization formalism of Sec. Ill Al 
We compose the scattering matrices Sl and Sr from the 
scattering matrix 



St — 



-^/l~T -iVT 



(48) 



of a semitransparent mirror with transmission probabil- 
ity T, the scattering matrix 

Sn. = ^— ( I'' ^ ) (49) 

' ni+n2 V2V^mHJ n2~ni J ^ ^ 

for a refractive index step from ni to n2, and the scat- 
tering matrix 

Sn^( . *^"P(*7^/2c)\ 

y exp{iu}nL/2c) J 



for ballistic propagation through a segment of refractive 
index n and length L/2. Then [s^ 



Sl — St ° •S'l^m, ° Sn^ ° Sn^,!-, 

Sr — Si^nn ° Snj^ O Sn^,! O St, 



(51a) 
(51b) 



with matrix composition rule ([T^ 

Following this prescription, we obtain the elements of 
5l in the form 



{riL + l)A+ + {UL - l)AlXl 

[nL + l)Al + {nT-l)Alxr ^ ^ 

^ 2inL\/TXL 

{nL + l)Al + {nL-l)Alxr 



{flL - l)Al + [UL + 1)A+X2 



{nL + l)Al + {nL-l)A+Xr 
where we introduced 

^1 



Al - -[(l + ani)Vl -T+l-dni], Xl ^ e 



Analogously, 



tr 



[UR - 1)A^ + {UR + 1)A+X^ 



t'R=tR = - 



{nR + l)A^ + {nR~l)A+X]^' 
2inR\/TXR 



{nR + l)A^ + {nR~l)A+Xl' 
{hr + 1)A+ + {hr - l)Aj^X^ 
{nR + l)A^ + {nR~l)A+X%' 



(52c) 

IL071L L/2c 

(53) 

(54a) 
(54b) 
(54c) 



where 



A]i^ -[{l + anR)VT^+l-cjnR], XR = e' 



jrL/2c 

(55) 

We now apply the quantization condition ([5]). This 
delivers the equation 



inL+nR)iAlA],-AtA],XiXf,) 
inL~nR){AlA+Xl-AlA],Xl), 



(56) 



which has to be solved for the resonance frequencies en- 
tering via Xl and Xr. 

Equation (j56p holds irrespective of whether the res- 
onator is T-'T-symmetric or not, but from now on we as- 
sume that this symmetry holds and therefore make use 
of ur — n\ = no — in J and A^ — [A^)* . Moreover, we 
introduce the scaled dimensionless frequency 

Q. — LonoL/c (57) 
and the dimensionless degree of nonhermiticity 

a — nj/nQ. (58) 
In terms of these quantities, we can rewrite Eq. (1561) as 



\At\ 



2 in 



\Al\^e~'^+ia{AlA+*e°'"^A+Al*e-°'^^) = 0. 

(59) 
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FIG. 2. Illustration of spontaneous PT-symmetry breaking 
in the model resonator of Fig. [ijc) , for T = (closed system) . 
The plot shows the dependence of the 12 lowest-lying dimen- 
sionless resonance frequencies fl on the hermiticity-breaking 
parameter a, following their trajectory as long as they are 
real. As a increases, adjacent resonance frequencies approach 
each other pairwise, until they merge in an exceptional point. 
Beyond the exceptional point, the involved frequencies be- 
come a complex-conjugated pair, and then are no longer plot- 
ted. These results are obtained by numerical solution of 
Eq. (|60l). 

1. Closed system 

In the limit T = of a closed system, = 1. The 

quantization condition (|59p then takes the form 

J'in) = a smh{aVL) + sin(f7) ^ 0, (60) 

which demands us to find the roots of a real function. 
The solutions are real, or occur in complex conjugated 
pairs, as required by VT symmetry. We denote these 
solutions as VLm (where m — 1,2, 3,...), so that we can 
refer back to them when we discuss the open system. 

The transition from a real to a complex spectrum is 
driven by the degree of nonhcrmiticity a. This transition 
involves exceptional points which occur when two real 
frequencies coalesce, delivering the additional condition 

^ — ^a^ cosh(al7) -I- cos{n) = 0. (61) 

For an exceptional point, Eqs. (|60)) and (|6T|) have to be 
fulfilled simultaneously. We denote the corresponding 
value of a and the value Vt of the two coalescing fre- 
quencies as al and Ct*, respectively, where Z = 1, 2, 3, . . . 
enumerates the exceptional points. Equations (1601) and 
([61]) then are equivalent to the conditions 

cos{VL*i) = -a*i, cosh{a*ini) = 1/af. (62) 

For illustration, we plot in Fig. [2] a set of quantized 
frequencies as a function of a, for a frequency range cov- 
ering the 12 lowest-lying levels. The figure displays six 



exceptional points at which pairs of real frequencies coa- 
lesce. Beyond these points, the involved frequencies be- 
come complex, and then are no longer plotted. 

For a = 0, the mth resonance frequency is located at 
f2m = mTT (m = 1, 2, 3, . . .). At the exception points two 
consecutive frequencies CL21-1 and Q,2i approach the value 
O* « 2Z7r — 7r/2, an expression which becomes more and 
more accurate as I increases. In the same limit, Eqs. ((60)) 
and (j6ip deliver the approximate condition 

Q,*garnr p. 2. (63) 

With increasing Z, decreases steadily, while the prod- 
uct increases very slowly. 

In order to investigate the behaviour close to an ex- 
ceptional point, we expand 

(64) 

where J- ~ dJ-/da, and we only kept the leading non- 
vanishing terms in this expansion. Slightly below the 
Ith exceptional point, the two resonance frequencies are 
therefore given by 



Because of the square-root dependence, the spectral ar- 
rangement near the exceptional point occurs over a very 
small range of a. In the following, we will use the term 
"far below (or above) the exceptional point" to indicate 
that the eigenvalues are well separated, while the regime 
where they are close together is called "slightly below 
(or above) the exceptional point" . However, in terms of 
numerical values a actually needs to approach a* very 
closely before the latter regime is entered. For later con- 
venience we also note that by a similar expansion as in 
Eq. ([M)) . the resonance splitting 

n2i{a) - n2i-i{a) « 2J-'(a, n2i)/:F"{ain^) (66) 

can be related to the deviation of from the condition 
(|6ip for an exceptional point. 

2. Slightly open system 

When the resonator is opened, bound states turn into 
quasibound states with complex resonance frequencies 
which are generally shifted by a mode- dependent amount 
downwards in the complex plane. For the system studied 
here, we find for T ^ 1 that this shift oc T is approx- 
imately rigid (i.e., mode independent), and is combined 
with a hfting - T^/^ of resonance coalescence close to 
the exceptional points (to obtain exact coalescence in 
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FIG. 3. Evolution of a pair of dimensionless resonance fre- 
quencies in the complex plane, while sweeping the nonher- 
miticity parameter a of the PT-symmetric model resonator 
depicted in Fig. [Tic). Shown are the resonance frequencies of 
the third and fourth level, with focus on the region around the 
exceptional point of the closed system (T = 0, as realized in 
panel a). Panels (6) and (c) show the lifting of exact degener- 
acy when the system is slightly open (T = 0.02 and T — 0.1, 
respectively). The arrows indicate the direction of the evo- 
lution with increasing a. Before the exceptional point, the 
resonances approach each other roughly parallel to the real 
frequency axis. Beyond the exceptional point, one of the two 
resonances is pushed towards the real axis, which facilitates 
the reaching of the lasing threshold (see Ref. [l3| for a system 
following the same scenario). These results are obtained by 
numerical solution of Eq. (f59|l, where we set no = 2. 



the complex plane, additional parameters beside a would 
have to be varied). These features follow by expanding 
the resonance condition (|59p up to order T^, upon which 
it takes the form 



where 



.To 



IYJ"(0) -a(l + a2) 



r2 

2^ °sinhar2 = 



To 



HqT 



l-T/2' 

To linear order in Tq, Eq. (|67p is solved by 

= — iro/2. 



(67) 



(68) 



(69) 



where r2,„ are the bound-state frequencies of the closed 
system, determined by Eq. (j60|) . Slightly opening up this 
resonator thus does not change the real parts of the res- 
onance frequencies, and shifts the imaginary parts by a 
mode- independent amount. In accordance to Eq. (|47l) . 
this shift does not depend on the degree of nonhermitic- 
ity a, which identifies Tq as the cold-cavity decay rate. 
Notably, Tq also determines the rigid shift of complex- 
conjugated pairs if T^T-symmetry in the closed system is 
spontaneously broken. 



For illustration, we plot in Fig. |3] the evolution of 
two resonance frequencies in the complex plane while 
the parameter a passes through an exceptional point 
of the closed system. The arrows indicate the evolu- 
tion direction of the resonance frequencies with increas- 
ing a. We set uq = 2, such that the expected shift 
ro/2 w r/(l - T/2) w T for smaU T. For the closed sys- 
tem (panel a), two real eigenvalues approach each other 
horizontally until they merge at the exceptional point, 
after with they become complex and move almost ver- 
tically away from each other. In panels (b) and (c), 
where the leakiness is small (T — 0.02) and moderate 
(T = 0.1), respectively, the resonance frequencies are 
shifted downwards by the expected amounts (for T = 0.1, 
r/(l - T/2) = 0.105, which corresponds well to the shift 
in the asymptotic regions to the left and right of the 
plotted range). The exceptional point is lifted only very 
slightly in panel (b), but much more distinctively in panel 
(c), and upon varying T we find that the distance of clos- 
est approach in the complex plane is indeed of order T'^/^. 



C. Frequency-resolved output radiation intensity 

We now turn to the frequency-resolved output radia- 
tion intensities Il{<-^) and Ir{lo) from the left and right 
openings of the resonator. These follow by substituting 
the refiection and transmission amplitudes (|52l) and ((54)) 
into Eqs. and (PSI) . respectively. Since these expres- 
sions contain the general resonance quantization condi- 
tion ([5]) in the denominator, the resonance frequencies 
calculated in the previous section determine the emission 
lines around which the output intensity is large. In this 
subsection, we directly derive expressions near resonance, 
including for the situation near an exceptional point; in 
the following subsection these are compared to the gen- 
eral near-resonant expression For the discussion, 
we again make use of the dimensionless variables fi and 
a [Eqs. (I57l) and (|58)) . respectively]. 

Before we present analytical results, we illustrate the 
key features in Fig. |4l where the solid lines are numerical 
results for I = II + obtained from the expressions 
([M)) and ([55]) as described above (the dashed lines rep- 
resent analytical results derived below). We set hq — 2, 
T — 0.02, and plot the intensity in a frequency range 
covering the third and the forth resonance, for values 
of a far below, near, and slightly beyond the excep- 
tional point 0-2 ~ 0.208573 of the closed system (reaching 
up very closely to the lasing threshold, which occurs at 
a « 0.208587). 

The upper panel in Fig. [4] shows results for four val- 
ues of a far below the exceptional point. The resonance 
frequencies are then well isolated, and the output inten- 
sity displays well-resolved resonance peaks which fit to 
a Lorentzian lineshape. As a increases, the resonances 
draw together and increase in height, but do not change 
their width, which is « Fq. The middle panel shows how 
resonances start to merge as one approaches the excep- 
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FIG. 4. (Color online) Frequency-resolved intensity = 
+ /^(i^) for the 'PT-symmetric model resonator de- 
picted in Fig. [TJc), for no = 2 and T — 0.02. Frequencies are 
in the range of the third and fourth resonance. Exact results 
from Eqs. (|34|) and Eq. ((35} (solid curves) are compared with 
various approximations (dashed curves), as the nonhermitic- 
ity parameter a is changed to steer the system from the regime 
of well-isolated resonances across the exceptional point of the 
closed system (q2 ~ 0.208573) up to just below the lasing 
threshold (the threshold itself is at a « 0.208587). Below the 
exceptional point (top panel), the isolated resonances agree 
with the Lorentzian lineshape ((71} . Slightly below the ex- 
ceptional point (middle panel), the interfering resonances are 
well described by Eq. ((75}. This expression also holds true 
for a = a2, where Eq. ((75} degenerates into the squared- 
Lorentzian (|73p (lower left panel). Close to the lasing thresh- 
old one obtains again a Lorentzian line shape, in agreement 
with (|76p (lower right panel). Note the dramatic increase in 
intensity as a is changed incrementally. 



tional point. This goes along with a dramatically increas- 
ing intensity, which is in accordance to the expectation 
from strongly-violated mode orthogonality. At the ex- 
ceptional point (lower left panel), the resonances merge 
into a very high single peak, which can be described by a 
squared Lorentzian, which still is of width Tq. Moving 



beyond the exceptional point, the resonance peak retains 
its centre. As one approaches the lasing threshold (where 
one of the complex resonance frequencies approaches the 
real axis) the emission line reverts to a Lorentzian, and 
the peak intensity increases further while the resonance 
width decreases. 

In order to explain these results we now obtain ana- 
lytical expressions for the resonance peaks in the output 
intensity, covering the whole range far below, near and 
slightly beyond the exceptional point (where one quickly 
reaches the lasing threshold). On obtaining these ex- 
pressions, it suffices to focus the attention on 7^, as Ir 
leads to the same results in the relevant leading orders 
in T. The key step in the derivation is to approximate 
the resonant denominator in Eq. (1341) by the same steps 
that lead from Eq. (j59p to the simplified quantization 
condition (|67p . Applying similar (more straightforward) 
approximations also to the denominator, this leads to the 
expression 



I Lin) 



1 En 



(1 + a^ysinhafl 



27r 2 



1^(0) 



a{l + a'^)-fsmhan\^ 
(70) 

from which all subsequent results follow. 

Eor an isolated resonance far below the exceptional 
point, we expand the denominator around the resonance 
condition fl = Clm, where is the real quantized fre- 
quency of the closed system, given by J" = [Eq. (jSU]) ]. 
This delivers the Lorentzian expression 



1 {l + a^ysmhan 



To/2 



2tt 



r^Ci^) (r!-an)2 + r2/4' 



(71) 

which in the upper panel of Fig. |4] is shown as a dashed 
line on top of the solid lines representing the exact numer- 
ical results. In the hermitian limit a — >■ 0, J^'{Clm) 1, 
and 



1 



aflmTo/2 



2tt (f7-an)2+r^/4 



(72) 



As expected, Eq. (|7T|) diverges at an exceptional point, 
which here is manifest because the term T'l^m) then 
vanishes [see Eq. ((6T|)]. This is remedied by keeping 
the next orders in the expansion of the denominator in 
Eq. (j70p . which results in the squared Lorentzian 



hm = 



1 {l + afysmha^n* 



2To 



2n 
1 



1 



2To 



271 af sinha^n;^ \ - fl^y + Tl/4\^ '' 



where we used Eq. (j60p to write 

F"{n*i)^a'i{l + af)sinhatni 



(73) 



(74) 



Equation (|73|) is shown in the lower left panel of Fig. [4] 
as the dashed curve on top of the numerical result at the 
exceptional point. 
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In order to better understand the approach to this 
squared Lorentzian, let us rederive it by considering the 
merging of the two involved resonance frequencies fl2i-i 
and CI21 of the closed system. While each associated par- 
tial intensity diverges as one approaches the exceptional 
point, their coherent sum 



1 To (l + af2)2sinhaff}f 



2n 2 



1 



1 



^ — ^21-1 
1 



- «ro/2 



n-^2i 
2ro 



zro/2 



(75) 

[where we have made use of Eqs. (1661) and (j74l) ] remains 
finite, and reduces to the squared Lorentzian ([75]) when 
the two resonance frequencies coalesce. Expression ([75)1 
is also accurate slightly away from the exceptional point, 
as can be seen from the curves in the middle panel of 
Fig.H 

Slightly beyond the exceptional point one approaches 
the lasing threshold, where one of the two resonances 
(with index 21 — 1 \i continuously labeled as in Fig. [3]) 
comes close to the real axis, thereby acquiring a much 
reduced width AJl = —2lTaVl2i-i ~ To — 2Imi72;-i <C 
r. The peak intensity of this resonance then exceeds by 
far that of the other resonance, resulting again in a (very 
narrow) Lorentzian 

J ^ 1 1 2/ro 

' 27raf sinhafr^f (r2-Rel72/_i)2 + Ar2V4' 

(76) 

This is plotted in the bottom right panel of Fig. |4l 



D. Implications of mode nonorthogonality 

We now discuss the preceding analytical results for the 
model resonator from the perspective of mode nonorthog- 
onality, based on the general considerations presented in 
Sec. HVl Our goal is to recover Eq. ([TTj) by substitut- 
ing the resonance wavefunction of the model system into 
Eq. (USD. 

Assuming no incoming radiation from the outside of 
the resonator, and employing scaled units s = x/L, the 
wave function in the left and right parts of the system 
(in which the respective refractive index is constant) can 
be written as 

V-L = r„^,Me'"(i+'")(^+i/2' + e-'f2(i+»a)(.+i/2)^ (77^) 
= a[e'"(i-^")(^-i/2) +^^^^_^^e-'^^(i-^°)(^-i/2)].(77b) 

Here 

rn^.M » -1 + nLT/2 - -1 + ^(1 + la), (78a) 
r„„,M«-l + nKr/2 = -l + ^^(l-za), (78b) 



are the reflection coefficients of the mirrors, including the 
step in the refractive index from 1 to or nn, respec- 
tively. The coefficient a is determined by the matching 
condition V'l(O) = V'fl(O), giving for Fq — > 



sinh[zf7(l + ia)/2] 
sinh[irj(l-m)/2]' 



(79) 



The second matching condition V'l(O) — '0i?(O) recovers 
the quantization condition (1601) . From this, we only need 
the previously established property f^™ w fi™ — iFo/2 
[Eq. (|M| ]. For Fq ^ and real 51^, we furthermore find 
|a| = l. 

Based on these expressions, and keeping only leading 
orders in Fq (which also appears in fim = ^l„i — iFo/2), 
as well as using the quantization condition ()60p . we find 



Imrj^n^l^Ap = -2n^(l 4- o?)Or,i sinh(aan), (80a) 



Imrj^n^lV-r = -2n^(l + «')f^mFo, 



7^+£ 



[2n2„n2^'(n,„)] 



(80b) 
(80c) 



Equation (|43p then delivers 



1 (H-a2)2sinh(ar2™) 



Fn 



27r 



j-'2(a„) (ri-i7„)2 + r2/4' 

(81) 

which indeed agrees exactly with our earlier result (j7ip . 
given that = Ir — 1/2 m the considered regime 
Fq ^ 1. It is noteworthy that we explicitly recov- 
ered the term J-' in the denominator, which vanishes 
at an exceptional point. This confirms that the inte- 
gral in the denominator of constitutes the appro- 
priate overlap integral, which quantifies the degree of 
mode nonorthogonality, as we already argued on general 
grounds in Sec. IIVI Finally, we once more call atten- 
tion to Eq. (j75p . which shows how the intensity is reg- 
ularized by the interference of the near-degenerate reso- 
nances. In this construction, the partial amplitudes are 
still consistent with Eq. (|5T|) . We thus can conclude that 
the appearance of the squared-Lorentzian line shape at 
an exceptional point can be explicitly linked to the self- 
orthogonality property of the resonance wave function, 
Eq. (Uni). 



VI. CONCLUSIONS 

Motivated by the recent interest in optical nonhermi- 
tian T'T-symmetric systems, we investigated the radia- 
tion intensity emitted by a resonator which is partially 
filled with an amplifying and an absorbing medium. This 
required to combine aspects of quantum noise with the 
properties of the resonator modes (especially, the con- 
sequences of broken hermiticity) . We addressed these 
issues on the common basis of scattering theory, which 
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allows to include quantum noise via the quantum-optical 
input-output formalism (see Sees. Ill Bl andllli p . while also 
giving access to the resonant frequencies (see Sec. Ill Ap 
and partial amplitudes in different regions of the res- 
onator (see Sec. IIVI) . 

Our main result is an expression of the near-resonant 
frequency-resolved radiation intensity, Eq. (j43|) . which re- 
lates this quantity to the properties of the resonant ra- 
diation mode, and includes an explicit measure of mode 
nonorthogonality induced by the amplifying and absorb- 
ing regions. T'T-symmetric systems provide natural ac- 
cess to exceptional points, where two resonances become 
degenerate, and their modes become self-orthogonal. The 
partial intensity of each resonance then diverges, but 
their sum yields a finite result, with a squared-Lorentzian 
line shape. Compared to the case of isolated resonances. 



we find that the total intensity is dramatically increased 
(see Fig. which should facilitate the observation of 
this radiation in experiments. 

We validated these results for the case of a model res- 
onator (see Sec. |V|, for which we obtained explicit an- 
alytical results in the whole physically accessible range 
of broken hermiticity, from the case of isolated reso- 
nances [Eq. ([7T|) ] over the near-degenerate case close to 
an exceptional point [Eq. (|75|) ] up to the lasing threshold 

[Eq. dzni)]. 
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